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Drawing the analogy to BaBiOa, we investigate via ab-initio electronic structure calculations 
potential new superconductors of the type ATIX3 with A = Rb, Cs and X = F, CI, and Br, with 
a particular emphasis on RbTlCla. Based on chemical reasoning, supported by the calculations, 
we show that Tl-based perovskites have structural and charge instabilities driven by the lone pair 
effect, similar to the case of BaBiOa, effectively becoming A2T1^"'"T1''"'"X6. We find that upon hole 
doping of RbTlCla, structures without Tl^^, Tl'^'^ charge disproportionation become more stable, 
although the ideal cubic perovskite, often viewed as the best host for superconductivity, should not 
be the most stable phase in the system. The known superconductor (Sr,K)Bi03 and hole doped 
RbTlCla, predicted to be most stable in the same tetragonal structure, display highly analogous 
calculated electronic band structures. 



I. INTRODUCTION 

Although not anticipated as such, relativistic effects 
play an important role in all parts of chemistry. Among 
the most prominent examples arc the color of An and 
the fact that Hg is a liquid at ambient conditionsii These 
observations are commonly attributed to the inert pair 
effect (IPE), which is the contraction and stabilization 
of the 6s orbitals in the context of the Lanthanide 
contraction and relativistic effects. The relativistic 
stabilization of the 6s orbitals greatly influences the 
chemistry of the heavy p-block elements. TICI3 (TP+ 
6s°), for example, readily decomposes to form TlCl (Tl^"'" 
6s^) and CI2. Relativistic effects are not only important 
in inorganic chemistry but also in materials science 
and condensed matter physics. The unique physics of 
topological insulators is due to compounds with heavy 
atoms that show strong relativistic effects<^ For a review 
on relativistic effects in general see reference 1. Here we 
study the relationship between superconductivity and 
charge instabilities due to the IPE in perovskite related 
structures. We focus on Tl compounds in particular, 
since Tl is the element with the most pronounced IPE 
known, and as a consequence is most stable in its +1 
(6s^) rather than its +3 (6s°) oxidation state. 

After the cuprates, the next highest Tc perovskite 
superconductor family is based on BaBiOs?^ which has 
a formal charge of +4 for Bi (which would be 6s^). 
From a chemical perspective, the correct formula for 
BaBiOs should be Ba2Bi2 06 with both and Bi5+ 

present. This 'charge disproportionation' is due to 
the IPE stabilization of the 6s^ configuration of Bi''^. 
BaBiOa crystallizes in a monoclinic distorted version of 
the perovskite structure*^ the BiOg-octahedra are tilted 
such that the angle between them is 159.9° in contrast to 
180° in the perfect perovskite structure. Due to the Bi'^+ 



and Bi'^"'' present, larger and smaller BiOg octahcdra 
are found. This distortion can also be described as 
an alternating breathing in and breathing out of the 
BiOg octahedra. The charge separation in BaBiOs is 
widely accepted in the litcraturc-S^— If BaBiOs is hole 
doped, either with K on the Ba site, or Pb on the Bi 
site, it becomes superconducting. The maximum Tc 
is 34 K for BaQ.6Ko.4Bi03 4 Bai^^K^BiOg (BKBO) 
crystallizes in the nominally simple cubic perovskite 
structure in the superconducting composition range 
0.375 < X < 0.5 without any evidence of a breathing 
mode distortion, suggesting that simple cubic symmetry 
is preferred for the superconductivityi^ A more recent 
report, however, claims tetragonal symmetry for BKBO, 
the distortion arising from rotations of the octahedraJ^ 
This is consistent with the structure of the Pb doped 
version. BaPbi-^BiajOa (BPBO) is superconducting for 
0.05 < X < 0.3 and its maximum Tc is 13 K for x = 0.25 
where there is a mixture of orthorhombic and tetragonal 
polymorphsfii In both systems, BKBO and BPBO, 
the maximum Tc is close to a metal to semiconductor 
transition. When BaBiOa is hole doped, more Bi^"'" 
is present and the average Bi-0 distance consequently 
decreases. Calculations have shown that a decreasing 
Bi-0 distance causes the band gap to decrease and 
leads to metallic conductivity when enough electrons 
are removedi^ This yields superconductivity at a critical 
concentration of Bi^"*". The idea that these types of 
valence instabilities could lead to superconductivity has 
long been arguedfi^iH and recent models quantify the 
presence of enhanced electron-phonon coupling due to 
this effect 



Here we show that the Tl perovskite RbTlCls, and by 
inference the larger class of hypothetical Tl perovskites 
(Rb,Cs)Tl(F,Cl,Br)3, display the same electronic fea- 
tures as BaBiOs. The nominal Tl^+ present is expected 
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to charge disproportionate into TP+ and Tl^+ due to 
the IPE, making RbTlCla a charge density wave (CDW) 
insulator hke BaBiOa and implying that superconduc- 
tivity should be possible with hole doping. The expected 
TcS of two hypothetical Cs variants with the ideal cu- 
bic perovskite structure in fact have recently been pre- 
dicted within the framework of the local density approx- 
imation (LDA), using linear response theory (LRT)ii^ 
Here we show that a van Hove singularity and a Fermi 
surface virtually identical to those seen in BaBiOa are 
present in the ideal cubic variant of the cubic Tl halide 
perovskites, but we also show by total energy calcula- 
tions that the ideal non-charge-disproportionated cubic 
structure should never be expected in these systems even 
when hole doped. Rather, lower symmetry perovskites 
with rotational distortions arc favored. Such phases are 
nonetheless also good hosts for superconductivity, as for 
example has been observed for (Sr,K)Bi03 

II. COMPUTATIONAL DETAILS 

The calculations were performed in the framework of 
density functional theory (DFT) using the wien2k— 
code with a full-potential linearized augmented plane- 
wave and local orbitals [FP-LAPW -|- lo] basisi^"— 
together with the Pcrdew Burke Ernzerhof (PBE) 
parameterization^! of the generalized gradient approx- 
imation (GGA) as the exchange-correlation functional. 
The plane wave cut-off parameter TImt^aiax was set 
to 7 and the irreducible Brillouin zone was sampled by 
40-100 k-points. For each structure type the lattice con- 
stants were estimated with the program SPuDSj^ and 
optimized by minimizing the total energy. The atomic 
positions were optimized by minimization of the forces. 

III. RESULTS 

The thallium halide perovskite CsTlCla is a known 
compound which was reported once in the literatures^ 
the other possible compounds are not so far reported. To 
estimate whether the formation of a perovskite should 
be possible, as well as its degree of distortion, the Gold- 
schmitt tolerance factor t is a good first step^^ where 
rx, TA and tm = anion radius, A site ion radius, and B 
site ion radius, respectively: 

t=-LA±I^ (1) 

V2[rM +rx] 

For a cubic structure, geometry requires that t = I. How- 
ever, the ideal cubic structure occurs for 0.89 < t < 1. 
Distorted perovskites, where the octahedra tilt and ro- 
tate around their shared corners to accommodate the size 
of the ion in the cavity, appear for 0.8 < t < 0.89. For 
smaller t's or t > 1, other structure types are usually fa- 
vorable. However, pervoskitcs are not truly ionic and the 
value of the tolerance factor depends on the value for the 



radii, which are not fixed for non ionic compounds. The 
tolerance factor can therefore only be used as a rough 
estimate. Using the Shannon radii,— we calculate the 
perovskite tolerance factor t (equation 1) for various Tl- 
based halide perovskites. The radius of Tl was estimated 
to be the average of the radius of Tl^+ and TP"*". The re- 
sults are shown in Table 1. For comparison, the tolerance 
factor for BaBiOa is 0.94. 

TABLE I. Tolerance factors for (Rb,Cs)Tl(F,CI,Br)3. 



Compound 


t 


RbTlFg 


0.86 


RbTlClg 


0.83 


RbTlBrs 


0.82 


CsTlFs 


0.90 


CsTlCls 


0.87 


CsTlBrs 


0.86 



The tolerance factors indicate that all of the com- 
pounds can crystallize in the perovskite structure, 
however, except for CsTlFs, they are all expected to 
be distorted through rotations of the octahedra. If 
the Tl compounds are hole doped on either the halide 
or the alkali metal site, the radius of Tl (rjvf) will 
become smaller because more TP"*" would be present. 
According to equation 1 the tolerance factor increases 
with decreasing rjvf . Hence higher symmetry structures 
should be possible with hole doping. 

We performed calculations for a model Tl halide per- 
ovskite, RbTlCla, in four different structure types: (a) an 
ideal undistorted cubic version (space group Pm3m, no. 
221), (b) a version with tilting of the octahedra as the 
only distortion (space group I4mcm, no. 140, the sym- 
metry of superconducting (Sr,K)Bi03i^) (c) a structure 
with breathing mode distortion without tilting (space 
groups Fm3m (225), the usual cubic double perovskite 
symmetry) and (d) one with both the breathing mode 
distortion and tilting of the octahedra (R3 (148), the 
symmetry of high temperature BaBiOg^) . All structures 
are presented in Figure 1. 

Comparison of the total energies of the four possi- 
ble RbTlCla structure types shows that structures with 
breathing mode distortions, i.e. the presence of IPE- 
driven charge disproportionation, arc much more stable 
than structures without breathing modes (Figure 2). As 
expected, we found such structures to be insulators with 
a band gap of the order of 1 cV. Structures without 
breathing modes calculate to be metallic without any 
doping, reinforcing the idea that the charge dispropor- 
tionation stabilizes the structure by opening a band gap 
(Figure 3). 

The band structure of the hypothetical undistorted 
ideal cubic structure of RbTlCla shows a distinct van 
Hove singularity (vHS) at the X point (see Figure 4a), 
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FIG. 1. The four different structure types of RbTlCls whose 
electronic structures are considered, a) the undistorted ideal 
cubic version, b) the tetragonal version with octahedra tilt- 
ing as the only distortion (no breathing mode distortion), c) 
the cubic version with breathing modes and no tilting, d) 
the rhombohedral structure with both distortions, tilting and 
breathing modes. Big spheres represent Rb atoms, medium 
one Tl atoms and small ones CI atoms. 




FIG. 2. Total energies of RbTlCls for different doping levels 
in the four structure types as a function of cell volume. 

which is highly analogous to the one seen in BaBiOa.— 
The presence of the vHS at Ej? is indicative of a struc- 
tural instability, where a distortion to decrease the den- 
sity of states at Ei? is strongly favored. The vHS has been 
cited as important for the occurrence of superconductiv- 
ity in perovskites and other materialsJ^ It also suggests 
that a structural distortion is necessary to stabilize this 
compound, driven by the need to decrease the density 
of states (DOS) at E^? that is present due to the vHS in 
the ideal cubic phase. The very strong electronic analogy 
of undoped ideal cubic RbTlCls to undoped ideal cubic 
BaBiOa can be seen by comparison of the Fermi surfaces 
in figures 4b and 4c. 

In order to simulate the effects of hole doping in 




E-Ep (eV) E-Ep (eV) 

FIG. 3. DOS of RbTlClg for different doping levels in the four 
structure types. 




FIG. 4. (a) Band structure and density of states of RbTlCls 
in the cubic undistorted structure. A distinct vHS is visible at 
the X point, (b) Fermi surfaces of RbTlCls and (c) BaBiOa. 



RbTlCla, we performed calculations within the virtual 
crystal approximation, with fewer electrons on the Rb 
site, allowing us to simulate the hypothetical hole-doped 
compounds Rbo.sTlCla and Rbo.eTlCla. We find that 
the volume decreases in each structure type in accor- 
dance with the increase of the Tl^+ concentration. Hole 
doping closes the band gaps in the semiconducting struc- 
tures. Comparison of the total energies shows that with 
increasing hole concentrations the energy differences be- 
tween the structures with and without breathing modes 
becomes smaller. Figure 2 shows plots of energy ver- 
sus volume for the different structures, for doped and 
undoped compounds. For the undoped compound, the 
breathing mode structures have a significantly lower en- 
ergy than structures with regular octahedra, even for the 
case of dramatically decreased cell volumes. Therefore 
an ideal cubic perovskite structure should not be sta- 
ble for undoped RbTlCls for experimentally accessible 
pressures, a reflection of the very strong IPE. With a 
hole doping level of 0.2 per formula unit, the tetragonal 
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structure with no breathing modes becomes lower in en- 
ergy, and, at hole dopings of 0.4 per formula unit, the 
tilted, non charge disproportionated structure is clearly 
the most stable. Therefore, with hole doping the breath- 
ing mode distortion becomes unfavorable, and a higher 
symmetry is preferred, although again the ideal cubic 
phase should not be seen. A similar scenario is observed 
in BKBO, where, at a doping level of 0.4 holes per for- 
mula unit, there is no evidence for breathing mode distor- 
tions, and a recent paper reveals a tilted-only structural 
distortion for the superconducting crystal structure In 
addition, we performed calculations for CsTlFs, the hy- 
pothetical perovskite in the thallium halide family with 
the highest tolerance factor (t = 0.9) to test whether 
it might adopt the ideal cubic structure with doping or 
pressure. However the energy curves were qualitatively 
similar to those seen for RbTlCla, from which we con- 
clude that the ideal cubic structure should not be stable 
for this family as a whole. 

In Figure 5 the DOS in the vicinity of for the 
four structures with and without doping is shown. The 
two structures with breathing mode distortions become 
metallic with doping, but the DOS at the Fermi level is 
very high. In the tetragonal structure, the DOS is much 
lower for Rbo.eTlCla, which again suggests the disappear- 
ance of breathing modes with doping. On doping, the 
energy gained by Tl disproportionation becomes smaller 
and the tilting of the octahedra becomes more impor- 
tant. This indicates that disproportionation is more fa- 
vorable if the average oxidation number is +2, i.e. that 
equal amounts of Tl^+ and Tl'^"'" are present and the ma- 
terial is insulating. Hole doping with x = 0.4 changes 
the average oxidation number to 2.4, which makes the 
stabilization by charge disproportionation less effective 
and should lead to identical bond lengths in the differ- 
ent TlClg octahedra. The increase of the unfavorable 6s° 
high oxidation state of TP+ could then lead to electronic 
instabilities that eventually induce superconductivity. 

To further draw the comparison between RbTlCla and 
the bismuthates, we calculated the band structure of 
Sro.4Ko.6Bi03, the composition with the highest Tc in 
the (Sr,K)Bi03 system, within the virtual crystal ap- 
proximation at the appropriate electron count Srg.yBiOa. 
Sro.4Ko.6Bi03 crystallizes in the same tetragonal crys- 
tal structure that we chose to model the tilting only 
structure for RbTlCls. Figure 5 shows the compari- 
son of the band structures of tetragonal Rbo.eTlCla and 
Sro.4Ko.6Bi03. The band structures are very similar, re- 
inforcing the strong electronic analogy between Bi-0 and 
Tl-Halide perovskites. 

IV. SUMMARY AND CONCLUSIONS 

Because the predicted most stable undopcd form of 
RbTlCla has frozen- in breathing modes, the scenario of 



charge disproportionation in Tl due to the IPE is realis- 
tic. The compound RbTlCla should be semiconducting if 

Rb.JlCL Sr„ RiO, 

0.6 3 0.7 3 
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FIG. 5. Comparison between the band structures of 
Rbo.eTlCls (left) and Sro.yBiOs (right). 



undoped, just like BaBiOa. Hole doping makes the band 
structure metallic and results in a stabilization of higher 
symmetry structures without frozen breathing modes, 
hence the possibility of superconductivity is present, in 
analogy to BaBiOa. The electron-phonon coupling, an 
additional important factor for the occurrence of super- 
conductivity, has also been predicted to be favorable in 
thallium halide perovskites^^, complementing this pic- 
ture. 

Other compounds in the family ATlXa with A = Rb, 
Cs, and X = F, CI, Br should have similar features. Fol- 
lowing the idea that the inert pair effect is the major 
reason for superconductivity in BKBO, one should also 
look for superconductivity in compounds in which Pb 
has an average oxidation number of +3. This could be 
present for example in CsPbCl20 or BaPb02F. 
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